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Recent preliminary Phenix data are consistent with a negative and sizable longitudinal double- 
spin asymmetry A^ h for 7T production at moderate transverse momentum p± ~ 1 4 4 GeV and 
central rapidity. By means of a systematic investigation of the relevant degrees of freedom we show 
that the perturbative QCD framework at leading power in p± produces at best a very small negative 
asymmetry in this kinematic range. 
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INTRODUCTION 



HARD-SCATTERING CALCULATION 



The determination of the nucleon's polarized gluon 
density is a major goal of current experiments with lon- 
gitudinally polarized protons at RHIC 1]. It can be ac- 
cessed through measurement of the spin asymmetries 



All = 



dAcr dcr++ - d(T+- 



da+- 



dcH 



(1) 



for high transverse momentum (p±) reactions such as 
pp — ► irX, jetX, jX. In Eq. QJ, <r ++ (cr + ~) denotes 
the cross section for scattering of two protons with same 
(opposite) helicities. High p^ implies large momentum 
transfer, and the cross sections for such reactions may 
be factorized into long-distance pieces that contain the 
desired information on the (spin) structure of the nu- 
cleon, and short-distance parts that describe the hard 
interactions of the partons and are amenable to QCD 
perturbation theory (pQCD). All of the reactions listed 
above have partonic Born cross sections involving gluons 
in the initial state and may therefore serve to examine 
the gluon content of the scattering longitudinally polar- 
ized protons. 

In this paper, we will consider the spin asymmetry 
A£ L for high-px tt° production, for which very recently 
the Phenix collaboration has presented first preliminary 
data at center-of-mass (c.m.s.) energy = 200 GeV 
and central rapidity. The data are consistent with a siz- 
able (up to a few per cent) negative A£ L in the region 
p± ~ 14-4 GeV. Even though the experimental uncertain- 
ties are still large and leave room for a different behavior 
of A£ L , trie new data motivate us to entertain the unex- 
pected possibility of A£ L being negative. To our knowl- 
edge, there have been no predictions of a substantially 
negative A^ L in the literature. This is in itself interest- 
ing, also because at moderate p± the spin asymmetry 
A£l i s expected to be particularly sensitive to gluon po- 
larization. 



We may write the polarized high-p^ cross section as 
/ dx a / dx b / dz c Aa(x a , fj,) Ab(xb, fi) 



dpj_dr] 



a,b,c 



dA ^ab(P-L,?7,3;a,Z&,Zc,M) 

dp±drj 



£>c(*cM) , (2) 



where r\ is the pion's pseudorapidity. The Aa, Ab (a, b = 
q, q, g) are the polarized parton densities; for instance, 



Ag(x,(j,) = g+{x,n) - g~{x,p) , 



(3) 



(the sign referring to the gluon helicity in a proton of 
positive helicity) is the polarized gluon distribution. The 
sum in Eq. (J5J) is over all partonic channels a + b — > c+ X, 
with their associated polarized cross sections dAa^ b . 
These start at 0{a 2 s ) in the strong coupling with the 
QCD tree-level scatterings: (i) gg — > gg, (ii) gg — > qq, 

(iii) gq(q) -> gq(q), ( iv ) qq -> ??, qq -> gg, qq -> qq, 

qq 1 — ► qq' , qq — > q'q' . The transition of parton c into 
the observed 7r° is described by the (spin-independent) 
fragmentation function D*. The functions in Eq. 
are tied together by their dependence on the factoriza- 
tion/renormalization scale \i which is of the order of the 
hard scale p±, but not further specified. All next- to- 
leading order [NLO, O(a^)] QCD contributions to po- 
larized parton scattering are known 3]. Corrections to 
Eq. J2J) itself are down by inverse powers of p± and are 
thus expected to become relevant if p± is not much bigger 
than typical hadronic mass scales. We neglect such con- 
tributions for now, but will briefly return to them later. 

To set the stage for our further considerations, Fig. 
shows NLO predictions for A£ L , for various gluon polar- 
izations Ag, all proposed within the framework of the 
analysis of data from polarized deeply-inelastic scatter- 
ing (DIS) in Q. Despite the fact that the Ag's used in 
Fig. Hare all very different from one another, none of the 
resulting A£ L is negative in the p± region we display. 
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FIG. 1: NLO predictions for A1 L based on different assum 
tions about Ag at the input scale for the evolution in 



BASIC OBSERVATIONS 

As a first step in our more detailed analysis of A£ L 
we discuss the various ingredients to Eq. J2J), other than 
the unknown Ag that one hopes to probe: the Ag and 
Ag distributions, the calculated partonic cross sections, 
and the fragmentation functions D™. They all jointly 
determine the analyzing power for Ag provided by A£ L . 

To a first approximation, the Ag, Ag may be treated 
as "known" . Ag + Ag has been probed extensively in 
polarized DIS. There certainly remains much room for 
improvement in our knowledge about them, but this does 
not affect A^ L too strongly. To give one example: for w° 
production the qg channel depends on the sum Au+Au+ 
Ad+Ad, so that uncertainties relating to the SU(2) flavor 
structure of the polarized sea are not relevant. 

We now turn to the partonic cross sections and frag- 
mentation functions. We focus on the lowest order (LO) 
which is expected to capture the qualitatively important 
points. Among the reactions (i)-(iii) listed above that 
have gluons in the initial state, process (ii) has a nega- 
tive partonic spin asymmetry oll = — lj while (i) and 
(iii) both have oll > 0. Also, for the p± region 
of interest, one has an average (z c ) ~ 0.3 -r- 0.4, where 
quarks are more favored than gluons to fragment into pi- 
ons. A first guess is, then, to attribute a negative A^ L 
to the negative gg — > gg cross section. However, this 
expectation is refuted by the numerical hierarchy in the 
partonic cross sections: at f\ = in the partonic c.m.s., 
which is most relevant for the Phenix data, channel (i) 
is (in absolute magnitude) larger than (ii) by a factor 
of about 160. Therefore, if one wanted to suppress the 
(positive) gg — > gg contribution, one would effectively 
have to switch off the gluon-to-pion fragmentation func- 
tion Dg. Even though our knowledge of Dg is incom- 
plete, this does not appear to be a sensible solution, for 
two reasons. First, data from e + e~ collisions, in par- 



ticular analyses of hadron production in bb plus (mostly 
gluon) jet final states []|, do constrain Dg significantly, 
even at fairly large z c . For example, the D^'s in the sets 
of |6| and [jj are in reasonable agreement with these data, 
with the one of arguably setting a lower bound on Dg . 
Second, elimination of the gg — > gg channel would also 
strongly affect the unpolarized cross section and reduce 
it by up to an order of magnitude at RHIC energies. 
However, previous Phenix measurements of the un- 
polarized cross section for pp — > n °A at VS = 200 GeV 
were found to be in excellent agreement with NLO cal- 
culations employing the of |6j . 

We therefore exclude that the gg — > qq channel is in- 
strumental in making v4£ L negative, and we thus have to 
investigate possibilities within Ag itself, and its involve- 
ment in gg — + gg and qg — > qg scattering. Given that the 
polarized scattering cross sections for these reactions are 
both positive, and that the first process comes roughly 
with the square of Ag, it is immediately clear that a siz- 
able negative asymmetry will not be easily obtained. In 
the next section, we will demonstrate this for the particu- 
larly instructive case of rapidity integrated cross sections. 



A LOWER BOUND ON A£ L 

We consider the LO cross section integrated over all 
rapidities r\. This is not immediately relevant for com- 
parison to experiment, but it does capture the main point 
we want to make. From Eq. J2Jl, one obtains: 



p^dAer^ 
dpj_ 



dx a Aa{x a ,pi) 
dz c D*(z c ,(j,) 



p 3 ± dAa 
dpj 



dxb Ab(xb,fi) 
« 6 (4,M),(4) 



where x\ = x?^/z 2 x a Xb, xt = 2p±/VS- It is then conve- 
nient to take Mellin moments in x\ of the cross section, 



Act* (TV) ee / dx 2 T (x 2 T ) 



2 i 2\N-i p 3 ± dAa* 



(5) 



(6) 



10 dpx 
One obtains (we suppress the scale \i from now on) 

Ao-^(N) = &a N+1 Aa'f D^ 2N + 3 , 

where the AavX are the i^-moments of the partonic 
cross sections and, as usual, f N = j^ ) dxx N ~ 1 f{x) for 
the parton distribution and fragmentation functions. We 
now rewrite Eq. JBJ in a form that makes the dependence 
on the moments Ag N explicit: 

A<j"(N) = (Ag N+1 ) 2 A N + 2Ag N+1 B N + C N . (7) 

Here, A N represents the contributions from gg — * gg and 
gg — ► gg, B N the ones from qg — > qg, and C N those from 
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the (anti)quark scatterings (iv) above; in each case, the 
appropriate combinations of Ag, Aq distributions and 
fragmentation functions are included. 

Being a quadratic form in Ag N+1 , Aa 7r (N) possesses 
an extremum, given by the condition 



A N Ag 



N+l 



-B 



N 



(8) 



We note in passing that the same equation may also be 
derived directly from Eq. (T4J by regarding the cross sec- 
tion as a functional of Ag, using a variational approach, 
and taking Mellin moments of the resulting stationarity 
condition. In the following we neglect the contribution 
from the gg — > qq channel which, as we discussed above, 
is much smaller than that from gg —> gg for the p±_ we 
are interested in. The coefficient A N is then positive, 
and Eq. (JSJ) describes a minimum of Aa^(N), with value 



(B N Y/A N + C 



(9) 



It is then straightforward to perform a numerical Mellin 
inversion of this minimal cross section: 



p 3 ± dAa v 



dp A 



(10) 



where T denotes a suitable contour in complex- N space. 

For the numerical evaluation we use the LO Aq, Aq of 
GRSV 0|, the Dl of 6], and a fixed scale n = 2.5 GeV. 
We find that the minimal asymmetry resulting from this 
exercise is negative indeed, but very small: in the range 
p± ~ 1 -7- 4 GeV its absolute value does not exceed 
1CP 3 . The Ag in Eq. JHJ that minimizes the asymme- 
try is shown in Fig. compared to Ag of the GRSV LO 
"standard" set |4|. One can see that it has a node and 
is generally much smaller than the GRSV one, except at 
large x. The node makes it possible to probe the two 
gluon densities in the 1717 term at values of x a , Xb where 
they have different sign, which helps in decreasing A£ L . 

Even though we have made some approximations in 
deriving the bound in Eq. we do believe that it 

exhibits the basic difficulty with a sizable negative ^4£ L 
at moderate p±: the fact that the cross section is a 
quadratic form in Ag effectively means that it is bounded 
from below. Note that this bound does not always imply 
that a negative A£ L is small: at higher p± it does allow 
a fairly large AJ L . 

In its details, the qualitative picture drawn by our 
study is of course subject to a number of corrections. 
First of all, we have integrated the cross section over all 
rapidities, whereas for the Phenix data \rj\ < 0.38. It is 
instructive to investigate the qualitative differences asso- 
ciated with this. Figure U| shows the polarized LO cross 
section at p± = 2.5 GeV versus the "distance" \x a — Xb\ 
in parton momentum fractions in Eq. @, for integra- 
tion over all r\ and for \r\\ < 0.38. Here we have used the 
GRSV polarized parton densities. The larger the rapidity 
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FIG. 2: Ag(x,n = 2.5 GeV) resulting from Eq. © (solid). 
The dashed line shows the GRSV LO "standard" Ag fl. 
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FIG. 3: dAcr^/dpx in bins of \x a - x b \. 

range probed, the more likely become collisions of par- 
tons with rather different momentum fractions. Indeed, 
the distribution for |?7| < 0.38 is steeper, implying that a 
node in Ag(x,n) will now be somewhat less efficient in 
promoting negative values for the asymmetry. 

In a more realistic calculation one would also prefer 
fi ~ p± to a fixed fi. Furthermore, since the Ag and Ag 
are coupled to Ag via evolution, any change of Ag will 
require a re-tuning of the Ag, Ag densities, so that the 
agreement with the polarized DIS data remains intact. 
Inclusion of the NLO corrections is important as well. 

All these points can be thoroughly addressed only in a 
"global" NLO analysis of the data, taking into account 
the results from polarized DIS as well. We will now report 
on such an analysis. Given that the data are still prelim- 
inary, this may seem premature. Our primary goal, how- 
ever, is to investigate whether the findings of our some- 
what idealized case, as summarized by Eqs. © and (fTTTfl . 
hold true in general. 
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"GLOBAL" NLO ANALYSIS 



CONCLUSIONS 



The main technical difficulty in a full global NLO anal- 
ysis of polarized DIS and pp — > ir°X data is the numerical 
complexity of the evaluation of the NLO corrections for 
the latter cross section. A convenient way to alleviate 
this problem was presented in [9j. Starting from Eq. J2J), 
one expresses the Aa, Ab by their Mellin inverses, e.g., 

Aa(x,p) = — [ dNx- N Aa N (u). (11) 



After interchange of integrations one obtains 
dAo-* 



dp±dr] 



= [ dN I dM Aa N (p)Ab M (p) 



a,b ,JLN 



(12) 



where the contain the partonic cross sections, the 
fragmentation functions, and all integrations over mo- 
mentum fractions, with the factors x~ and x b as 
complex "dummy" parton distributions according to 
Eq. (|11|) . The strength of this approach is that there 
is no dependence of p* b on the moments Aa N , Ab M of 
the true parton densities. This means that the p* b can 
be pre-calculated prior to the analysis on a specific array 
of the two Mellin variables N and M. One chooses a 
convenient functional form for the parton distributions, 
depending on a set of free parameters. The latter are 
then determined from a x 2 minimization procedure. The 
double inverse Mellin transformation which finally links 
the parton distributions with the pre-calculated p^ b of 
course still needs to be performed in each step of the fit- 
ting procedure, but becomes extremely fast by choosing 
the values for N,M on the contours IV, Tm simply as 
the supports for a Gaussian integration. 

Following these lines, we have performed a simultane- 
ous analysis of all data from polarized DIS and of the 
preliminary Phenix data for pp — > tt°X. We have used 
several different functional forms for the polarized gluon 
density, in particular allowing it to have a node. The 
quark densities were allowed to vary as well. We have 
artificially decreased the error bars of the data points 
for A£ L in order to see whether the fit can be forced to 
reproduce a negative of about —2% in the region 
p± ~ 1 -r 4 GeV. We have also slightly shifted individual 
data points to study the response of the fit. In no case 
have we been able to find a fit that yielded a negative 
with absolute value larger than a few times 10~ 3 . 
Even those fits, however, gave a negative A£ L only at 
the higher end of the p± interval, and invariably they 
led to a polarized gluon density that had a node and 
tended to violate positivity \Ag\ < g in certain ranges of 
x. The "global" analysis thus confirms our qualitative 
finding above that any negative A^ L is also very small. 



Our analysis demonstrates that pQCD at leading 
power in p± predicts that is bounded from below 
by AJl ~ £>(-l(T 3 ) in the region p± ~ 1 4- 4 GeV. 
The observation relies on collinear factorization and on 
exploring the physically acceptable ranges of parton dis- 
tribution and fragmentation functions. 

For now, the data |3| do not allow a compelling con- 
clusion on whether the bound is violated or not. What 
should one conclude if future, more precise, data will 
indeed confirm a sizable negative A£ L ? As indicated ear- 
lier, corrections to Eq. J2J are associated with power- 
suppressed contributions to the cross section. Since p± 
is not too large, such contributions might well be sig- 
nificant. On the other hand, comparisons of unpolar- 
ized 7T° spectra measured at colliders with NLO QCD 
calculations do not exhibit any compelling trace of non- 
leading power effects even down to fairly lowpj_ > 1 GeV, 
within the uncertainties of the calculation. It is con- 
ceivable that the spin-dependent cross section with its 
fairly tedious cancelations has larger power-suppressed 
contributions than the unpolarized one. One may at- 
tempt to model the effects by implementing an "intrin- 
sic" transverse- momentum (fcj_) smearing for the initial 
partons which generically leads to corrections by powers 
of (k±)/p±, with (k±) an average k±. Such effects were 
shown to have indeed some potential impact on at 
p± < 5 GeV ^3 ■ A negative A£ L would open up a quite 
unexpected window on aspects of nucleon structure and 
limitations of pQCD so far little explored. 
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